Viral population estimation using pyrosequencing

The diversity of virus populations within single infected hosts presents a major difficulty for the natural immune response as well as for vaccine design and antiviral drug therapy. Recently developed pyrophosphate based sequencing technologies (pyrosequencing) can be used for quantifying this diversity by ultra-deep sequencing of virus samples. We present computational methods for the analysis of such sequence data and apply these techniques to pyrosequencing data obtained from HIV populations within patients harboring drug resistant virus strains. Our main result is the estimation of the population structure of the sample from the pyrosequencing reads. This inference is based on a statistical approach to error correction, followed by a combinatorial algorithm for constructing a minimal set of haplotypes that explain the data. Using this set of explaining haplotypes, we apply a statistical model to infer the frequencies of the haplotypes in the population via an EM algorithm. We demonstrate that pyrosequencing reads allow for effective population reconstruction by extensive simulations and by comparison to 165 sequences obtained directly from clonal sequencing of four independent, diverse HIV populations. Thus, pyrosequencing can be used for cost-effective estimation of the structure of virus populations, promising new insights into viral evolutionary dynamics and disease control strategies.Comment: 23 pages, 13 figure

Measuring Enzymatic HIV-1 Susceptibility to Two Reverse Transcriptase Inhibitors as a Rapid and Simple Approach to HIV-1 Drug-Resistance Testing

Simple and cost-effective approaches for HIV drug-resistance testing are highly desirable for managing increasingly expanding HIV-1 infected populations who initiate antiretroviral therapy (ART), particularly in resource-limited settings. Non-nucleoside reverse trancriptase inhibitor (NNRTI)-based regimens with an NRTI backbone containing lamivudine (3TC) or emtricitabine (FTC) are preferred first ART regimens. Failure with these drug combinations typically involves the selection of NNRTI- and/or 3TC/FTC- resistant viruses. Therefore, the availability of simple assays to measure both types of drug resistance is critical. We have developed a high throughput screening test for assessing enzymatic resistance of the HIV-1 RT in plasma to 3TC/FTC and NNRTIs. The test uses the sensitive “Amp-RT” assay with a newly-developed real-time PCR format to screen biochemically for drug resistance in single reactions containing either 3TC-triphosphate (3TC-TP) or nevirapine (NVP). Assay cut-offs were defined based on testing a large panel of subtype B and non-subtype B clinical samples with known genotypic profiles. Enzymatic 3TC resistance correlated well with the presence of M184I/V, and reduced NVP susceptibility was strongly associated with the presence of K103N, Y181C/I, Y188L, and G190A/Q. The sensitivity and specificity for detecting resistance were 97.0% and 96.0% in samples with M184V, and 97.4% and 96.2% for samples with NNRTI mutations, respectively. We further demonstrate the utility of an HIV capture method in plasma by using magnetic beads coated with CD44 antibody that eliminates the need for ultracentifugation. Thus our results support the use of this simple approach for distinguishing WT from NNRTI- or 3TC/FTC-resistant viruses in clinical samples. This enzymatic testing is subtype-independent and can assist in the clinical management of diverse populations particularly in resource-limited settings

Cross-validated stepwise regression for identification of novel non-nucleoside reverse transcriptase inhibitor resistance associated mutations

<p>Abstract</p> <p>Background</p> <p>Linear regression models are used to quantitatively predict drug resistance, the phenotype, from the HIV-1 viral genotype. As new antiretroviral drugs become available, new resistance pathways emerge and the number of resistance associated mutations continues to increase. To accurately identify which drug options are left, the main goal of the modeling has been to maximize predictivity and not interpretability. However, we originally selected linear regression as the preferred method for its transparency as opposed to other techniques such as neural networks. Here, we apply a method to lower the complexity of these phenotype prediction models using a 3-fold cross-validated selection of mutations.</p> <p>Results</p> <p>Compared to standard stepwise regression we were able to reduce the number of mutations in the reverse transcriptase (RT) inhibitor models as well as the number of interaction terms accounting for synergistic and antagonistic effects. This reduction in complexity was most significant for the non-nucleoside reverse transcriptase inhibitor (NNRTI) models, while maintaining prediction accuracy and retaining virtually all known resistance associated mutations as first order terms in the models. Furthermore, for etravirine (ETR) a better performance was seen on two years of unseen data. By analyzing the phenotype prediction models we identified a list of forty novel NNRTI mutations, putatively associated with resistance. The resistance association of novel variants at known NNRTI resistance positions: 100, 101, 181, 190, 221 and of mutations at positions not previously linked with NNRTI resistance: 102, 139, 219, 241, 376 and 382 was confirmed by phenotyping site-directed mutants.</p> <p>Conclusions</p> <p>We successfully identified and validated novel NNRTI resistance associated mutations by developing parsimonious resistance prediction models in which repeated cross-validation within the stepwise regression was applied. Our model selection technique is computationally feasible for large data sets and provides an approach to the continued identification of resistance-causing mutations.</p

Density‐ and size‐dependent mortality in fish early life stages

The importance of survival and growth variations early in life for population dynamics depends on the degrees of compensatory density dependence and size dependence in survival at later life stages. Quantifying density‐ and size‐dependent mortality at different juvenile stages is therefore important to understand and potentially predict the recruitment to the population. We applied a statistical state‐space modelling approach to analyse time series of abundance and mean body size of larval and juvenile fish. The focus was to identify the importance of abundance and body size for growth and survival through successive larval and juvenile age intervals, and to quantify how the dynamics propagate through the early life to influence recruitment. We thus identified both relevant ages and mechanisms (i.e. density dependence and size dependence in survival and growth) linking recruitment variability to early life dynamics. The analysis was conducted on six economically and ecologically important fish populations from cold temperate and sub‐arctic marine ecosystems. Our results underscore the importance of size for survival early in life. The comparative analysis suggests that size‐dependent mortality and density‐dependent growth frequently occur at a transition from pelagic to demersal habitats, which may be linked to competition for suitable habitat. The generality of this hypothesis warrants testing in future research.publishedVersio

From sea monsters to charismatic megafauna: changes in perception and use of large marine animals

Marine megafauna has always elicited contrasting feelings. In the past, large marine animals were often depicted as fantastic mythological creatures and dangerous monsters, while also arousing human curiosity. Marine megafauna has been a valuable resource to exploit, leading to the collapse of populations and local extinctions. In addition, some species have been perceived as competitors of fishers for marine resources and were often actively culled. Since the 1970s, there has been a change in the perception and use of megafauna. The growth of marine tourism, increasingly oriented towards the observation of wildlife, has driven a shift from extractive to non-extractive use, supporting the conservation of at least some species of marine megafauna. In this paper, we review and compare the changes in the perception and use of three megafaunal groups, cetaceans, elasmobranchs and groupers, with a special focus on European cultures. We highlight the main drivers and the timing of these changes, compare different taxonomic groups and species, and highlight the implications for management and conservation. One of the main drivers of the shift in perception, shared by all the three groups of megafauna, has been a general increase in curiosity towards wildlife, stimulated inter alia by documentaries (from the early 1970s onwards), and also promoted by easy access to scuba diving. At the same time, environmental campaigns have been developed to raise public awareness regarding marine wildlife, especially cetaceans, a process greatly facilitated by the rise of Internet and the World Wide Web. Currently, all the three groups (cetaceans, elasmobranchs and groupers) may represent valuable resources for ecotourism. Strikingly, the economic value of live specimens may exceed their value for human consumption. A further change in perception involving all the three groups is related to a growing understanding and appreciation of their key ecological role. The shift from extractive to non-extractive use has the potential for promoting species conservation and local economic growth. However, the change in use may not benefit the original stakeholders (e.g. fishers or whalers) and there may therefore be a case for providing compensation for disadvantaged stakeholders. Moreover, it is increasingly clear that even non-extractive use may have a negative impact on marine megafauna, therefore regulations are needed.SFRH/BPD/102494/2014, UID/MAR/04292/2019, IS1403info:eu-repo/semantics/publishedVersio

Integrated Moloney murine leukemia virus DNA studied by using complementary DNA which does not recognize endogenous related sequences.

By preannealing a radioactive, representative Moloney murine leukemia virus (M-MuLV) cDNA with large excesses of AKR 70S viral RNA, an M-MuLV-specific cDNA has been prepared. When hybridized to restriction enzyme fragments of M-MuLV-infected mouse cell DNA, the preannealed probe recognizes integrated M-MuLV DNA and does not recognize endogenous related DNA sequences found in uninfected mouse cells. The viral DNA sequences recognized by the preannealed probe are spread throughout the viral genome, although some sequences are recognized less efficiently. By using this preannealed probe, multiple integrations of M-MuLV DNA have been detected in infected fibroblasts and in an M-MuLV-induced tumor. Integrated viral DNA fragments smaller than the complete viral genome have also been detected. By using this preannealed probe to examine a mass-infected culture of mouse fibroblasts, no evidence for a strongly preferred site for M-MuLV integration could be found